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Amendment to the Claims : 

This listing of claims replaces all prior versions, and 
listings, of claims in the application: 

1. (Original) A method for estimating the frequency of a 
single frequency complex exponential tone in additive Gaussian 
noise, comprising the steps of: 

performing the fast Fourier transform (FFT) on the tone; 

estimating the frequency as the frequency corresponding to 
the largest FFT output coefficient magnitude; 

computing a discriminant which is proportional to the 
frequency error in the initial frequency estimate using modified 
coefficients of the discrete Fourier transform (DFT) with center 
frequencies plus one half and minus one half of the FFT bin 
spacing relative to the initial frequency estimate; 

mapping the value of the discriminant into the estimate of 
the frequency error in the initial frequency estimate using a 
mathematically derived function; 

adding the estimate of the frequency error to the initial 
frequency estimate to get a first interpolated frequency 
estimate ; 

computing a further discriminant which is proportional to 
the frequency error in the first interpolated frequency estimate 
using modified coefficients of the discrete Fourier transform 
(DFT) with center frequencies plus one half and minus one half 
of the FFT bin spacing relative to the first interpolated 
frequency estimate ; 

mapping the value of the further discriminant into the 
estimate of the frequency error in the first interpolated 
frequency estimate using the mathematically derived function; 
and 

adding the estimate of the frequency error in the first 
interpolated frequency estimate to the first interpolated 
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frequency estimate to get a second interpolated frequency 
estimate . 

2. (Original) The method according to claim 1, wherein the 
first interpolated frequency estimate is in a region of 
relatively low noise induced frequency error. 

3. (Currently amended) The method according to claim [1 or] 
2, wherein the method is implemented in computer hardware and/or 
computer software. 

4. (Currently amended) The method according to [any one of 
the preceding] claim [s] 1, wherein the method is utilised in 
communications, signal processing and biomedical applications. 

5. (Currently amended) The method according to [any one of 
the preceding] claim [s] 1, further comprising the steps of: 

iteratively deriving an interpolated frequency estimate, and 
using the frequency discriminant, to obtain a more precise 
frequency estimate. 

6. (Original) The method according to claim 5, wherein the 
steps of iteratively deriving an interpolated frequency estimate 
and using the frequency discriminant are repeated until a fixed 
point solution occurs, where at this fixed point, the 
discriminant function has zero value. 

7. (Currently amended) The method according to [any one of 
the preceding] claim [s] 1, wherein the frequency discriminant is 
computed by: 

D(M)='P'-' a ' 
\fi\ + \a\ 
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where, £ = fT s — ,rax 



N ' 

* = fT s -^, and 

□ and □ are the modified DFT coefficients defined by, 
P = Ydc™ + h = Xr(n)e ,2OT ^ 

Z n=0 

and, cc = Y(k max --) = Y i r(n)e i2m N 

^ n=0 

thus, the initial frequency estimate using the FFT, 

f 0 Tc = ^-and£ = 0. 
0 s N 

8. (Currently amended) The method according to [any one of] 
claim [s] 1 [to 6], wherein the frequency discriminant is 
computed by: 

D .li/»r-i«r >fory>0 ., 

Y\PY+\ccY 

where □ and □ are the modified DFT coefficients defined by, 
1 n^i _- 2m 2_ 



y9 = Y(k max+ -) = X r ( n ) e 



n=0 



vi i (^max ~) 

1 ^1 -j2 2_ 



and , a = - -) = £ r(n)e 

^ n=0 

9. (Currently amended) The method according to [any one of] 
claim [s ]1 [to 6], wherein the discriminant of frequency 
estimation error is computed by: 

i i/?i 2 -i<*i 2 

2|/?| 2 +|a| 2 ' 

where □ and □ are the modified DFT coefficients defined by, 
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1 £4 2_ 



y» = Y(k„ BX +-) = X r (n)e 



n=0 

and, a = YCk^ - ^) = j>(n)e 

^ n=0 

10. (Currently amended) The method according to [any one 
of] claim [s] 1 [to 6], wherein the frequency discriminant is 
computed by: 

D.Ra^] 
/3 + a 

where Re [ . ] is the real part and * denotes the complex 
conjugate, and □ and □ are the modified DFT coefficients defined 
by, 



/? = Y(k niax +-) = Xr(n)e 



n=0 



1 N^l _ j2/m 2_ 



and , a = YCk^ --) = £r(n)e 



n=0 



11. (Currently amended) The method according to [any one of 
preceding] claim [s] 1, wherein the frequency discriminant is 
computed by using more than two DFT coefficients. 

12. (Currently amended) The method according to claim 11, 

N 

wherein 2M+2 coefficients are used, where 0<M< 1 and the FFT 

2 

coefficients are used in the frequency discriminant with optimal 
weighting coefficients obtained by using the concept of matched 
filtering is, 
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Ml 1 

£C , {YtOc^ +- + m)modN]- Y'tOc^ ---m)modN]} 

D =Re[ r ^ \ 1 

£{C • < Y K k ™x +- + m)modN] + Y*[(k nm ---m)modN]} 

i^ + r m] 2 2 

N 

where, 0 <M < 1, mod N indicates modulo N, 

. 2 

and, where, * denotes complex conjugate. 

j/r[— +m][ ] 

_e 2 N sin|>(| + m)] 
2 sin[— (j + m)] 

and, YCk^ +-^ + m)and Y(k niax -^--m) are the modified DFT coefficients 
given by, 



Y(k max +i + m) = j>(n)e 



n=0 



1 ^1 _ j2 2 



and # Y(k m - - - m) = £ r(n)e 



n=0 



13. (Currently amended) The method according to claim [11 
or] 12, wherein the frequency discriminant using more than two 
DFT coefficients [us] d.s used in the last iteration to obtain 
additional frequency accuracy. 

14. (Currently amended) The method according to [any one 
of] claim[s] 11 [to 13], wherein the frequency discriminant is 
computed by using more than two DFT coefficients and less or 
equal to all N FFT coefficients. 

15. (Currently amended) The method according to [any one of 
the preceding] claim [s] 1, wherein additional frequency accuracy 
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is obtained by computing the frequency discriminant recursively 
until convergence for the frequency estimate is reached. 

16. (Original) The method according to claim 15, wherein 
convergence for the frequency estimate is reached after zero to 
three iterations, the number of iterations being dependent on 
the specific discriminant used and the signal to noise ratio. 



17. (Currently amended) The method according to claim [15 
or] 16, wherein in any iteration, the frequency discriminant is 
computed using any one of the functional forms: 

Af m (r) = ltan- 1 [ l/?ml " |aml tan(^)] f s , or 

Af m (r) =-±— [ { f m ~ | a ™ I" ] f s , where y m is a constant,^ > 0 , or 
Af m (r) = — [^iziEmi]f s t f or y=1/ or 



2N |/? m | + |a m | 
4N |AJ 2 +|a n 



*m(r)~l \^i m 2 ]f s / for y=2 



18. (Original) The method according to claim 17, wherein □ 
varies on each iteration. 

19. (Currently amended) The method according to claim [15 
or] 16, wherein in any iteration, the frequency discriminant is 
computed using: 

Afm(r) = ^r Re[ * m ]f * ' 

2N p m +a m 

where, Re [ . ] denotes the real part and * denotes the complex 
conjugate . 



20. (Currently amended) The method according to [any one 
of] claim [s] 17 [to 19], wherein the frequency incremental 
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shift, Af m (r) , is related to the previously defined frequency 
discriminant, D, by, 

Af m (r) = -^D 
m ' 2N 

21. (Currently amended) The method according to [any one 
of] claimts] 15 [to 20], wherein the frequency discriminant is 
driven to zero input and output values by either modifying the 
frequency of the DFT coefficients or frequency translating the 
signal . 

22. (Currently amended) The method according to [any one 
of] claimts] 15 [to 21], wherein signal frequency translation is 
achieved, by multiplication of the signal by a locally generated 
complex exponential signal. 

23. (Original) The method according to claim 22, wherein 
frequency multiplication of the signal is implemented with a 
standard hardware, software, or combination hardware/software 
FFT . 

24. (Original) The method according to claim 23, wherein the 
hardware/software FFT is highly optimized for at least one 
processor operating as a system. 

25. (Currently amended) The method according to [any one 
of] claimts] 15 [to 24], further comprising the step of scaling 
the frequency estimate during recursion, to save multiplies. 

26. (Original) The method according to claim 25, further 
comprising a final step of multiplying the scaled frequency 
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estimate. f m+ iT s with the sampling frequency f s to remove the 
scaling from the frequency estimate. 

27. (Currently amended) A frequency estimation software 
program for estimating the frequency of a single frequency 
complex exponential tone in additive Gaussian noise, wherein the 
frequency estimation program has functionality to perform the 
method according to [any one of] claim [s] 1 [to 26] . 

28. (Currently amended) A computer system programmed to 
perform the method according to [any one of] claim [s] 1 [to 26] . 

29. (Currently amended) The computer system according to 
claim 28, wherein the hardware includes a DSP processor chip. 
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